Seeing zeros of random polynomials: quantized vortices in the ideal Bose gas 
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We propose a physical system allowing one to experimentally observe the distribution of the 
complex zeros of a random polynomial. We consider a degenerate, rotating, quasi-ideal atomic 
Bose gas prepared in the lowest Landau level. Thermal fluctuations provide the randomness of the 
bosonic fleld and of the locations of the vortex cores. These vortices can be mapped to zeros of 
random polynomials, and observed in the density proflle of the gas. 
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An important field of study in theoretical statistical 
physics concerns the properties of the roots of random 
polynomials P, Q] ■ Of particular interest is the so-called 
Weyl polynomial for the complex variable C,: 



Z^n vm! 



(1) 



where the a^'s are independent random complex num- 
bers with the same Gaussian probability distribution. 
The roots of P in the complex plane can be mapped 
to a two-dimensional (2D) gas of particles with repulsive 
interactions. They are spatially antibunched, and have a 
uniform mean density in the large Wmax limit . 

Although the statistical properties of the roots of the 
Weyl polynomial have been well studied theoretically, no 
physical system has allowed yet to observe them directly. 
The goal of this Letter is to show that a 2D rotating 
ideal Bose gas is a well suited system for this observa- 
tion. The positions of the vortices appearing in the gas 
can be mapped to the zeroes of the random polynomial 
describing the atomic state. More precisely, the gas is 
harmonically trapped and it rotates at a frequency close 
to the trapping frequency so that it is 'frozen' in the 
lowest Landau level (LLL) 0, 13- Its finite temperature 
T ensures that several vortices are present in the region 
where atomic density is significant, and thermal fluctua- 
tions provide the randomness of the vortex locations in 
different realizations of the experiment. Such an experi- 
ment is not unrealistic: a 2D atomic Bose gas in the LLL 
has recently been produced 0, and the use of a Fano- 
Feshbach scattering resonance allows to nearly cancel the 
interactions between ultra-cold atoms 0] ■ This regime of 
an ideal gas with large thermal fluctuations dramatically 
differs from the well studied case of a rotating and in- 
teracting Bose-Einstein condensate in the LLL at T = 0, 
where the vortices are known to form an Abrikosov lat- 
tice 111113. 

The ideal gas in our model is confined in a harmonic 
trap, with oscillation frequency cu in the xy plane and 
ujz along z. The confinement along z is assumed to be 
strong, k^T ^ ^icj^, so that the z degree of freedom is 



frozen and the gas is kinematically 2D. We also assume 
that some angular momentum has been transferred to 
the gas by a stirring procedure , so that thermal equi- 
librium is reached in a frame rotating at frequency 
around z. is chosen close to the trapping frequency lu: 



u; — <C a; 
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which is experimentally realistic since the value (to — 
il)/uj ~ 0.01 has already been achieved 0. We also as- 
sume that the gas is cooled to a low enough temperature, 
fceT ^ 2fiuj, so that the relevant single particle_states are 
linear combinations of the LLL eigenmodes: 



C = X + iy 



(3) 



Here aho = ^Jh/Mu) (where M is the atomic mass) is 
taken as the unit of length. The mode energy = 
mh{uj — n) depends on a single quantum number m > 
so that thermodynamically the gas is effectively ID. 

The relevant quantity in our study is the density of 
vortices, which we will define in relation with the complex 
classical field ^(r) describing the state of the gas. This 
classical field represents not only the condensate (atoms 
in the vortex-free mode 0o) but also all superimposed 
thermal fluctuations. We obtain ip using the expression 
for the many-body density operator a of the ideal gas in 
thermal equilibrium in the grand canonical ensemble 0| : 
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In this expression ct is a statistical mixture of Glauber 
coherent states |coh : with positive weights -P({'0}) 
(the so-called Glauber P distribution) given by the Gaus- 
sian functional specified below. A given realization of the 
experiment can then be viewed as a random draw of the 
atomic field state ip, with the probability law P{{iJ;}). 

The stochastic nature of ip is simple to characterize by 
expanding it on the eigenmodes (j)rn' 



V'(r) 



m>0 
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The flm's are complex, statistically independent random 
numbers, with a Gaussian law: 
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where Um = [exp(/3(em — /i)) — 1]^^ is the mean occu- 
pation number of mode m. Here /? = and is 
the chemical potential. This provides a clear link with 
the random polynomial of Eq. |^ when several rim have 
similar values: 



V'(r)=/(C)- 
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When /(C) is factorized as /(C) cx ni(C~Ci)j each root Q 
corresponds to the location of a positively charged vortex 
in the field "0. Since having a multiple root is a zero 
measure event, these vortices are of charge unity. Note 
that ip results from the interference of a large number 
of macroscopically populated field modes, reminiscent of 
the interference of independent condensates 0, ^| . 

The standard case where all the Om's have the same 
variance corresponds in our model to all em's being equal, 
i.e. — uj. The average vortex density is then uni- 
form, py — 1/tt and the average pair distribution function 
P2{r — u/2, r -I- u/2), which depends only on the relative 
distance u, can be calculated analytically 0- However 
this case can not be achieved experimentally since for 
^ = CO, the centrifugal force exactly balances the trap- 
ping force and the gas is not confined anymore. In a 
realistic model one must address the case fl < uj, for 
which the trapping force overcomes the centrifugal one. 
The statistics of the roots then do not coincide with the 
standard results of the literature, and we must perform 
a study of this specific model. 

To provide an intuitive understanding on how the roots 
of /(C) arc distributed, we first show in figure ^numer- 
ical results for a randomly generated field ip. We take 
kgT = 500 ft(cij — f2), which is compatible with the con- 
dition k^T <C 2hLu for the experimentally realistic value 
U = 0.999 cj. On the first and second lines of figure^ we 
show plots of Indi/ip) and \tp\'^, respectively, for 3 values 
of p. In all cases, the roots locations are clearly visible. 

We now turn to an analytic study of the problem and 
we express the algebraic density of vortices pv (r) in terms 
of ip and its derivatives ■ A helpful simplification is 
to eliminate the Gaussian factor in Eq.(|7I), which does 
not change the vortex locations, and use /(C) rather than 
'!/)(r) as a random field. Using dyf{() = idxfiC) = «/'(C)i 
we obtain: 



p.(r) = |/f 5(2)[/(C)] 



(8) 



making it obvious that all vortices have a positive charge 
in the LLL. 

The expectation value of p^ gives the average vortex 
density py, the correlation function (pt,(r — u/2)p„(r -I- 
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FIG. 1: Density plots for three randomly generated fields 
■i/), using the grand canonical ensemble for ksT = 500 ei 
(ei = h{u3 — fl)). Left column: p = — lOei, mean number 
of particles A'^ = 1986; middle column: p = — ei, A'' — 3396; 
right column: p = -0.1 ei, = 8320. First line: ln(|V'(r)p). 
Second line: |i/;(r)p. Third line: simulation of a real exper- 
iment; random atomic positions are generated according 
to the density |'(/)(r)p, and a 129x129 pixel 'camera' image is 
produced assuming that each atom produces a Gaussian spot 
of (7 = 0.2aho. The width of each image is lOflho- Experimen- 
tally, a time of flight technique can be used to magnify these 
images, since the LLL modes are self-similar under a ballistic 
expansion. 



u/2)) gives the vortex pair distribution function p2, etc. 
Using properties of Gaussian statistics, these quantities 
can be expressed in terms of expectation values of prod- 
ucts of two fields, such as /, /*, or their derivatives with 
respect to C The exphcit result for p^ is: 



Pv{r) 
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The expectation values (. . .) are readily calculated using 
{a^Um') — rimSm.m'- Figurc [21 givcs the variation of p„ 



with r, for the temperature and the three values of the 
chemical potential used in figure ^ 

The general expression for the pair distribution func- 
tion p2 is rather involved, and we give it only for small 
relative distances it. It vanishes quadratically with m, 
showing the effective repulsion between the vortices: 



P2(r- 



u 



r+-^) =C(r)^ + 0(u3) 



(10) 



where C is a function of /, /*, and their first and second 
derivatives [l^ . 
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We find that the experimentally relevant situation fl < 
uj leads to results similar to those of the standard random 
polynomial theory (formally corresponding to = lu) in 
the vicinity of the trap center when the two following 
conditions are fulfilled. Firstly, many eigenmodes have 
to be thermally populated: 



(11) 



Secondly, the low energy modes must have compara- 
ble populations, imposing a negative chemical potential 
^ h{uj — n). This is the situation depicted in the left 
column of figure ^ and in figure |2^. 

If \fj,\ is reduced for a fixed T, the mean number of 
atoms increases and an imbalance appears among the 
populations of the low energy modes. The vortex density 
is then depleted near the trap center, as it can be seen 
in the middle and right columns of figure ^ and in figure 
|2b,c. Eventually, for ^ h{uj — ^l), the total population 
of the modes m 7^ saturates and a condensate forms in 
the mode m ~ 0. This condensate has no vortex and it 
expels the thermal vortices, which accumulate in a corona 
with a density locally exceeding 1/tt (figure . 

The situation at the onset of condensation is relevant 
experimentally since it is tempting to increase the num- 
ber of atoms in order to improve the experimental sig- 
nal. To extend quantitatively our analysis to this regime, 
we must replace the grand canonical ensemble, which 
presents non physical, large fluctuations of the total par- 
ticle number when a condensate is present, by the canoni- 
cal ensemble. We then adapt our definition of the vortex 
distribution by projecting Eq.Q) on a subspace with a 
fixed total number of particles N, resulting in the follow- 
ing probability distribution of the field tp: 



2N 



m 



(12) 



where = / IV'P- The key change with respect to 

the grand canonical ensemble is that the probability dis- 
tribution of the field is no longer Gaussian, showing that 
the proposed physical system opens a new class of prob- 
lems. We evaluate numerically the vortex density using 
a generating function technique In the trap center, 
one can also construct an exact mapping to a solvable 
problem, the calculation of the canonical partition func- 
tion and occupation numbers for a ID ideal Bose gas in 
a harmonic trap . This leads to the exact expression: 



PS(r = 0) = 



1 - e^/^^i 



e/3€i _ 1 



(13) 



where nf is the mean number of particles in the first 
excited mode m = 1 (with ei = h{uj — fi)) in the 
canonical ensemble As shown in figure the 

correct canonical result drops to a much smaller value 




FIG. 2: Vortex (p„) and atomic (p) densities. Solid line: 
result for pv in the grand canonical ensemble. Dasfied line: 
result for p„ in the canonical ensemble with the same mean 
number of particles A'^. Dotted line: result for p/100 in the 
grand canonical ensemble. fcsT — 500 ei with ei = h{uj — Q). 
(a) II = -10 ei {N = 1986). (b) /i = -ei (iV = 3396). (c) 
fj, = -0.1 ei (iV = 8320). The unit length is aho- 



than the incorrect grand canonical result in the trap cen- 
ter. In the large N limit, one can show that pl{0) ~ 

exp[— (A^ -I- l)/3ei]/7r whereas the grand canonical pre- 
diction tends to zero only as 1/iV. Note that p^(0) drops 
with N similarly to the probability exp(— 7V/3ei) of hav- 
ing an empty condensate mode, which is a natural con- 
dition to have a vortex in r = 0. 

How to observe the vortex density in practice? In cur- 
rent experiments with rotating interacting condensates, 
one measures the positions of the particles and the vor- 
tices appear as holes in the density profile . This strat- 
egy can be used in regions of space where the mean den- 
sity of particles p greatly exceeds the density of vortices 
fjy, here ~ I/tt. To each vortex embedded in such a high 
density region will correspond a clearly identifiable hole 
in the particle distribution. The maximal density of the 
non-condensed fraction of the gas is ^ kBT/[TT h{uj — il)] 
close to the trap center, obtained for <^ h{uj — fl). 
Then p ^ py when Ea. (|ll|l is satisfied. To indicate up to 
what distance from the trap center this condition holds, 
we have also plotted p(r)/100 in figure [3 

Conversely, do all the holes embedded in high density 
regions correspond to vortices? They may in principle 
correspond to a local minimum of where the field 
assumes a small but non zero value. Such 'spurious' local 
minima can in general form in a non-condensed ideal Bose 
gas, which is subject to large density fiuctuations due to 
the thermal bunching effect of the bosons. However, we 
can show that this phenomenon is absent in the LLL. 
We Taylor-expand the field to second order in u = 
r — ro around the location Tq of a stationary point of 
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l^/'P, assuming that ilj{ro) ^ 0: 

|^(r)|2 = |V;(ro)P [1 + u • Afu + 0(u=^)] 

where the 2x2 matrix M is real symmetric, 
finds that the trace of M is 



Tr M = Re 



grad 



(14) 
One then 



(15) 



where we used Eq.JT)) and grad(|?/'p) = to get the last 
identity. Since its trace is < 0, M cannot be positive and 
1^1^ cannot have a local minimum at Vq. 

We exemplify the above discussion by a Monte Carlo 
simulation of a real experiment for the parameters of fig- 
ure ^ The results are shown on the last line of figure ^ 
Starting from the random fields tp previously generated, 
we have produced these images by generating random po- 
sitions of particles according to the distribution and 
by mimicking the finite imaging resolution of a real exper- 
iment. Local minima of the density are visible, and can 
be checked on the images of the first and second lines to 
correspond to vortices. The visibility of the vortex pat- 
tern could be improved by increasing further the atom 
density. For a given ratio k-QT/hw, this can be achieved 
by rotating the gas even faster, so that (a; — il)/tj de- 
creases. 

In real life, there are interactions between the particles, 
characterized by the 3D s-wave scattering length a. The 
interaction potential in the LLL is then modeled by the 
pseudo-potential g5'^'^\Yi — r2), with g = l^ph^hwaj 
where az = \JhjMjjJz. We now derive a condition on a 
for the interactions to be negligible, in calculating to first 
order in g the effect of the interactions on the probabil- 
ity distribution of the vortices. Focusing on the quasi- 
homogeneous regime \n\ ^ hiio — fi), we set f2 = w. 
The unperturbed occupation numbers then all have the 
same value riQ. The vortex density being uniform, even 
in presence of interactions, we consider the first order 
correction in g to the vortex pair distribution function 
P2(u). Its expression can be simplified in the degenerate 
regime tt-q ^ 1: 



'?P2(u) 
pf(u) 



4 / d\ 



(p,(0)p.(u)|^(r)|^ 
pf(u) 



(ivxr)r 



(16) 

where the average is taken over the unperturbed distri- 
bution. The right hand side of Eq. (|16|l can be expressed 
analytically. It is maximal (in absolute value) in u = 0: 



lim 

M— »0 



(^102 (u) 

pf(u) 



IGtt 



9P_ 



(17) 



whereas it drops as u^eT^ at large distances. Interac- 
tions will then play a negligible role in the vortex distri- 
bution if the mean field energy is much smaller than 
1^1 . For the realistic numbers ksT = hio/2, uq — 100 



and ttz = 1 /zm, the condition for negligible interactions is 
\a\ < 0.2 nm. The scattering length can indeed be tuned 
in practice to such a low value using a Fano-Feshbach 
resonance. This tuning should be done only during the 
last fraction of the experimental sequence, after having 
taken advantage of a larger value of a to reach thermal 
equilibrium in a reasonable time. 

To summarize, the rotating quasi-ideal Bose gas is a 
promising system to implement in practice the concepts 
developed in the theory of random polynomials. It also 
raises novel questions such as the influence of the non- 
Gaussian statistics for the polynomial coefficients when a 
condensate is present. Finally we emphasize that our dis- 
ordered vortex pattern appears in a non-superfluid sys- 
tem. An interesting future line of research is the tran- 
sition to an ordered Abrikosov lattice when interactions 
are increased and the system becomes superfluid. 
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